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METHOD OF AND APPARATUS FOR MEASURING 
THICKNESS OF THIN FILM OR THIN LAYER 

The present invention relates to a method of and apparatus for measuring the 
thickness of a thin film or thin layer using the spectroscopic measurement. The method and 
apparatus can be used in various fields including, for example, the measurement of the 
thickness of various layers on a wafer substrate in the semiconductor industry. 

BACKGROUND OF THE INVENTION 

An application of a spectrophotometer using visible light or infrared light is the 
measurement of the thickness of a thin film or thin layer on a substrate. The principle of the 
spectroscopic thickness measurement is as follows. 

When a ray of incident light (or measurement light) 10 having a single wavelength is 
cast on a thin layer (or sample) S, as shown in Fig. 8, a part of the light is reflected on the 
front surface SI of the sample S, and the remainder enters the sample S. A part of the 
incoming light is also reflected by the rear surface S2 (i.e. the boundary with the substrate), 
goes back through the sample S, and goes out through the front surface SI. Since the first 
reflected light Rl and the second reflected light R2 have different optical path lengths, an 
interference occurs between the two light waves Rl and R2, depending on the wavelength X 
of the measurement light 10 and the thickness d of the sample S. When a graph is drawn 
with the wavelength (or wavenumber) of the measurement light as the abscissa and the 
intensity of the interference light as the ordinate while the wavelength of the measurement 
light is changed (or scanned), a wavy interference spectrum is obtained. The waveform of 
the interference spectrum can be represented by a cosine function whose cycle interval 
corresponds to the thickness d of the sample. Therefore, using the interference spectrum, it is 



possible to determine the thickness d of the sample S by the following steps: automatically or 
manually measuring the wavenumber at each crest (or peak) or trough (or valley) of the 
interference spectrum, determining the cycle interval in the wavenumbers between the crests 
or between the troughs by the least square error method or some other methods, and 
calculating the thickness d from the wavenumber cycle interval and a known refractive index 
n. 

The interference spectrum obtained through the spectroscopy rarely takes the ideal 
shape due to various factors, such as the wavenumber dependency of the interference 
efficiency, the wavenumber dependency of the energy distribution of the light source used, 
and various noises arising from the apparatus. Conventional methods do not take such 
factors into account; the methods assume that the waveform of the interference spectrum 
becomes an ideal cosine curve. Thus, it has been difficult to improve the accuracy of the 
thickness measurement. 

In view of the above problem, the applicant has proposed a method of measuring the 
thickness of a thin film or thin layer in the Japanese Patent Application No. 2002-147107. 
For a spectrum obtained by a measurement and represented by a graph with the wavenumber 
of the incident light as the abscissa and the intensity of the interference light as the ordinate, 
the method defines an approximate spectrum, called the "constructed spectrum," by 
modifying the ideal cosine curve while taking into account various factors that disturb the 
waveform. The constructed spectrum is represented by a function with thickness as a 
variable. Then, a graph is created to show the least square errors between the measured 
spectrum and the constructed spectrums for various thicknesses, and the thickness 
corresponding to the minimum point of the least square error is determined as the desired 
thickness. By taking into account the factors that disturb the waveform of the interference 
spectrum, the above method has improved the accuracy of the thickness measurement. 



In the process of determining the minimum point of the least square error between 
the constructed spectrum and the measured spectrum by the method of the aforementioned 
Japanese Patent Application, the thickness d is changed as a parameter to give various 
waveforms to the constructed spectrum. This method is applicable to the measurement of the 
thickness of not only a single-layer film but also a multiple layered film. 

The measurement of a multiple layered film, however, is possible only when all the 
layers have the same refractive index n. 

SUMMARY OF THE INVENTION 

The present invention addresses the problem, and an object of the present invention 
is, therefore, to provide a method of and an apparatus for measuring the thickness of a thin 
layer, which can measure the thickness of each layer of a multiple layered film even when 
the layers have different refractive indices. 

Thus, the present invention provides a method of measuring the thickness of every 
layer of a thin film consisting of a single or multiple layers, which include the steps of: 

casting a measurement light onto the thin film; 

observing a measured spectrum of an interference light of lights reflected by a front 
surface of the thin film, reflected by each boundary between the layers of the thin film, and 
reflected by a rear surface of the thin film, wherein the measured spectrum shows a cyclic 
curve; 

creating a constructed spectrum including a parameter, wherein the parameter has the 
same dimension as that representing the cyclic curve of the measured spectrum; 

calculating the least square error between the measured spectrum and the constructed 
spectrum; 

detecting a local minimum value of the least square error while changing the 



parameter within a predetermined range; and 

calculating the thickness of each layer of the thin film from a parameter yielding the 
local minimum value, using refractive indices of substances constituting the layers. 

In the above method, the measured spectrum can be observed using a wavelength 
dispersing element, such as a diffraction grating or a prism, and a linear array of photodiodes 
arranged along the dispersion of the wavelength. 

The present invention also provides an apparatus for measuring a thickness of every 
layer of a thin film consisting of a single or multiple layers, which includes: 

means for casting a measurement light onto the thin film; 

means for observing a measured spectrum of an interference light of lights reflected 
by a front surface of the thin film, reflected by each boundary between the layers of the thin 
film, and reflected by a rear surface of the thin film, wherein the measured spectrum shows a 
cyclic curve; 

means for creating a constructed spectrum including a parameter, wherein the 
parameter has the same dimension as that representing the cyclic curve of the measured 
spectrum; 

means for calculating the least square error between the measured spectrum and the 
constructed spectrum; 

means for detecting the local minimum value of the least square error while changing 
the parameter within a predetermined range; and 

means for calculating the thickness of each layer of the thin film from the parameter 
yielding the local minimum value, using refractive indices of the substances constituting the 
layers. 

Another method of the present invention includes the steps of: 
casting a monochromatic measurement light onto the thin film; 



observing a measured spectrum of an interference light of lights reflected by a front 
surface of the thin film, reflected by each boundary between the layers of the thin film, and 
reflected by a rear surface of the thin film while changing a wavelength of the measurement 
light, wherein the measured spectrum shows a cyclic curve; 

creating a constructed spectrum including a parameter, wherein the parameter has the 
same dimension as that representing the cyclic curve of the measured spectrum, 

calculating the least square error between the measured spectrum and the constructed 
spectrum, 

detecting a local minimum value of the least square error while changing the 
parameter within a predetermined range, and 

calculating the thickness of each layer of the thin film from the parameter yielding 
the local minimum value, using refractive indices of substances constituting the layers. 

Another apparatus of the present invention corresponding to the above method is an 
apparatus including: 

means for casting a monochromatic measurement light onto the thin; 

means for measuring an intensity of an interference light of lights reflected by a front 
surface of the thin film, reflected by each boundary of the layers of the thin film, and 
reflected by a rear surface of the thin film; 

means for creating a measured spectrum with the measured intensity of the 
interferance light while changing the wavelength of the measurement light; 

means for creating a constructed spectrum including a parameter, wherein the 
parameter has the same dimension as that representing the cyclic curve of the measured 
spectrum, 

means for calculating the least square error between the measured spectrum and the 
constructed spectrum, 



means for detecting a local minimum value of the least square error while changing 
the parameter within a predetermined range, and 

means for calculating the thickness of each layer of the thin film from the parameter 
yielding the local minimum value, using refractive indices of the substances constituting the 
5 layers. 

In the method or apparatus according to the present invention, when the sample film 
is a multiple layered film, the constructed spectrum is created by a linear sum of plural base 
spectrums, each having a cycle interval as a parameter, where the dimension of the cycle 
interval is the same as that of the cycle interval of the measured spectrum. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 shows an example of a multiple layered film onto which a measurement light 
is cast and by which the light is reflected. 

Fig. 2 shows an example of the measured spectrum of a two-layered film. 
1 5 Fig. 3 shows the relation between the measured vector and the constructed vector. 

Fig. 4 shows an example of a two-layered film onto which a measurement light is 
cast and by which the light is reflected. 

Fig. 5 shows the relation between the cycle interval and the least square error 
between the measured spectrum and the constructed spectrum obtained by the measurement 
20 of a two-layered film. 

Fig. 6 schematically shows the construction of an apparatus embodying the present 
invention. 

Fig. 7 shows the steps of calculating the thickness of a layer by the apparatus 
embodying the present invention. 
25 Fig. 8 illustrates the principle of the spectroscopic measurement of the thickness of a 
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film. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
A method embodying the present invention is described. 
5 Suppose a measurement light having a single wavelength is cast onto a multiple 

layered sample film 10 composed of i layers having the refractive indices of m, n 2 , ... , n i9 
respectively, as shown in Fig.l. Then, the measurement light is reflected by the front surface 
(i.e. the boundary surface between the film and the air), the boundary surfaces between the 
layers, and the rear surface (the boundary surface between the film and the ambient air) of 

10 the sample film 10. Thus, the sample film 10 produces (z+1) rays of reflected light with the 
(/+1) layers of the boundary surfaces. The reflected light waves interfere with each other, 
producing an interference light. The interference light is detected as the reflected light from 
the surface of the sample film 10. 

Changing (or scanning) the wavelength of the measurement light will change the 

1 5 intensity of the reflected light coming from the sample film 10. When a graph is drawn with 
the wavenumber or wavelength as the abscissa and the intensity of the reflected light as the 
ordinate, the graph becomes a cyclic spectrum, which is hereby called the "measured 
spectrum." Fig. 2 shows an example of the measured spectrum obtained by the 
measurement of a film composed of two layers having different refractive indices and being 

2 0 200nm and 4000nm in thickness. 

The following description assumes that the wavenumber is taken as the abscissa. In 
this case, the measured spectrum will show the following characteristics (or factors). 

(i) There is a certain regular cycle. 

(ii) There is an offset. 

2 5 (iii) There is an almost linear increasing drift. 



(iv) The amplitude of the interference wave decreases as the wavenumber increases 
because the interference efficiency decreases as the wavenumber increases. 

The characteristic (i) probably results from the interference of two light waves. A 
film having i layers provides (i+l) reflective surfaces. In this case, the number of different 
5 combinations of two interfering light waves is (,+i)C2. Denoting ( /+ i)C 2 by j, the characteristic 
(i) probably results from the superimposition of the j pieces of cyclic interference spectrums 
each produced by a different pair of light waves. For example, a film having two layers 
produces a reflected light composed of three cyclic interference spectrums superimposed on 
each other. The characteristics (ii) are probably due to the background noise from the 
1 0 apparatus. 

The measured spectrum shown in Fig. 2 seems to be composed of a first cosine wave 
having the cycle interval of about 1000 cm -1 , a second cosine wave having the cycle interval 
more than ten times as long as that of the first cosine wave, an offset and an increasing drift, 
all superimposed on each other. The offset and the drift are commonly observed components, 

15 as explained above. Since the sample film in Fig. 2 is composed of two layers, there must be 
three cyclic spectrums superimposed on each other. However, there are only two cyclic 
spectrums observable in the graph of Fig. 2. This fact shows that the thickness of each layer 
of the film cannot be directly calculated from the measured spectrum. In such a case, it is 
necessary to use the method according to the present invention. 

20 Letting g>i, a>2 ? CDj denote the angular frequencies of the two-beam interference 

spectrums of the aforementioned characteristic (i), and taking the other factors (ii), (iii) and 
(iv) into consideration, the spectrum can be approximately represented by the following 
equation: 

f(x)=a 0 +a i x+a 2 ( 1 /x) sin(co i x+8 1 )^ 

25 ...(1), 
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where the angular frequencies a>i, . . ., C0j have the dimension of length because, in the present 
case, the abscissa of the measurement spectrum represents the wavenumber, which is 
denoted by "x." 

In the right side of equation (1), the first term represents the offset, the second term 
5 represents the drift, and the third and subsequent terms represent cyclic waveforms. More 
precisely, the (1/x) of each of the third and subsequent terms shows the decrease in the 
amplitude according to the increase in the wavenumber, and 5 means the shift in the phase, 
which becomes conspicuous as the thickness becomes very large. 

Since 

1 0 sin(cox + 8) = sincaxcosS + coscoxsin5, 

equation (1) can be rewritten as 

f(x) = ao + (Xix + a2(l/x)sin(coix) + 013(1 /x)cos(coix) + ... + a2j(l/x)sin(a)jx) + 
a2j+ 1 ( 1 /x)cos(cDjx) . . . (2) . 

Using (2j+2) functions f 0 (x), fi(x), . . ., f 2j (x) and f2j+i(x) defined by 
15 fo(x) = l, 

fi(x) = x, 

f 2 (x) = (l/x)sin(coix), 
f 3 (x) = (l/x)cos(coix), 

2 0 f 2 j(x) = ( 1 /x)sin(cojx), and 

f 2 j+i(x) = (l/x)cos((Ojx), 
the function f(x) can be written as a linear sum of the above (2j+2) functions, as follows: 

f(x) = aofo(x) + a t fi(x) + a 2 f 2 (x) + a 3 f3(x) + . . . + a 2 jf 2 j(x) + a 2 j+if2j+i(x) ... (3). 

As explained above, the function f(x) is an appropriate approximation of a measured 
25 spectrum. Therefore, the following description uses the function f(x) as the constructed 
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spectrum corresponding to a measured spectrum g(x). The function f(x) has j unknown 
values coi, coj. 

The square error between the function f(x) and the measured spectrum g(x) will be 
minimized when the function f(x) is the best approximation of the measured spectrum g(x). 
Taking this into account, a set of values is given to the unknown angular frequencies [coi,co 2 , 
©j], and the coefficients a 0 , ai, a 2j+ i are changed so that the square error between the 
function f(x) and the measured spectrum g(x) is minimized. This operation is performed for 
various sets of angular frequencies [coi,a> 2 , • ^jL ^ d the least square error is obtained for 
each set of angular frequencies [coi,co 2 , coj]. The angular frequencies [a>i, co 2 , coj] 
that give the minimum value of the least square error is chosen as the most appropriate 
approximation of the angular frequencies of the j cyclic waveforms superimposed on each 
other in the measured spectrum. 

The above calculation, however, is impractical because it deals with j unknown 
values [coi, co 2 , coj]. Taking this into account, it is assumed here that there is only one 
unknown angular frequency, co. Using four functions fao(x), fai(x), fa 2 (x) and fa 3 (x) defined 
by 

fao(x)= 1, 
fai(x) = x, 

fa 2 (x) = (l/x)sin(cox), and 
fa 3 (x) = (l/x)cos(o>x), 

the function fa(x) is written here as a linear sum of the above four functions, as follows: 

fa(x) - aofao(x) + ajfai(x) + a 2 fa 2 (x) + a 3 fa 3 (x) ... (3'). 
The equation (3') has only one unknown frequency value, o. With the frequency co fixed at a 
certain value, the square error between the function fa(x) and the measured spectrum g(x) is 
calculated while changing the coefficients a 0 , a u a 2 and a 3 . The function fa(x) that gives the 
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minimum value of the square error is selected as the constructed spectrum for the given 
frequency co. Performing this operation for various co values, it is possible to find a value of 
go that minimizes the least square error. The frequency co thus found corresponds to the 
function fa(x) that is the best approximation of the measured spectrum g(x). 
5 Now, suppose the error between the measured spectrum and the constructed 

spectrum is represented by an error function whose variables include the angular frequency 
co and the cycle interval T. Here, the measured spectrum and the constructed spectrum are 
each represented by a vector whose elements include the intensity corresponding to 
wavenumbers predetermined at intervals of AT, or corresponding to wavelengths 

1 0 predetermined at intervals of AL 

As explained above, the constructed spectrum fa(x) is defined here as a linear sum of 
the four functions fao(x), fai(x), fa 2 (x) and fa 3 (x). This means the four functions fao(x), fai(x), 
fa 2 (x) and fa 3 (x) can be regarded as four base vectors whose vectorial sum provides the 
constructed vector fa(x). In this case, the distance between the measured vector Y and the 

15 subspace W defined by the four base vectors fao(x), fai(x), fa 2 (x) and fa 3 (x) in the vector 
space corresponds to the local minimum value of the least square error between the 
measured spectrum g(x) and the constructed spectrum fa(x) (see Fig. 3). That is, the best 
constructed spectrum fa(x) of the measured spectrum corresponds to the constructed vector 
Y' of the foot of the perpendicular line dropped from the measured vector Y onto the 

2 0 subspace W. The constructed vector Y' can be calculated as follows. 

The measured vector Y and the vector X representing the measurement point 
(wavenumber) for the measured vector Y are given as follows. 
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With the four base vectors fao(x), fai(x), fa2(x) and fa3(x) expressed using the elements of the 
vector X, a matrix A(T) is defined as the function of the cycle interval T=2nl(S), as shown by 
the following equation (4). 
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As explained above, the constructed vector Y' is defined by the foot of the 
perpendicular line dropped from the measured vector Y onto the subspace W. Therefore, the 
vector Y'-Y is orthogonal to each of the four base vectors. This relation can be expressed by 
the equation t A(Y'-Y)=0, which can be rewritten into the following equation (5) to obtain 
the constructed vector Y\ 

Y'^A^AaY'AY ...(5) 

Using the measured vector Y and the constructed vector Y\ the error function e(T) is 
defined here by the following equation. 



sin-^x N _, N-l 



...(4) 
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...(6) 



Equation (6) provides the least square error between the measured spectrum and the 
constructed spectrum for a given cycle interval T. When a graph is drawn with the cycle 
interval T as the abscissa and the error function e(T) as the ordinate, the graph will have one 
or more local minimum points. Each cycle interval T at which each local minimum point 
appears corresponds to one of the cycle intervals of the j modes of two-beam interferences. 

The j modes of two-beam interferences can be divided into two groups. The first 
group includes i modes of interferences, each of which is produced by two-beams coming 
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from the front and rear surfaces of each of the i layers. The other group includes (/-/) modes 
of interferences, each of which is produced by two-beams coming from the front and rear 
surfaces of two or more layers stacked. For example, the interference between the reflected 
light waves Rll and R12 in Fig. 1 requires only the first layer. Similarly, the interference 
between the reflected light waves R12 and R13 requires only the second layer. The 
interference between the reflected light waves R12 and R13, on the other hand, requires the 
first and second layers. The thickness of each of the / layers can be obtained by analyzing 
only the aforementioned i kinds of interferences. 

Thus, in practical measurements, it is necessary to identify the local minimum points 
that correspond to the interferences each produced by each of the / layers. This is practically 
possible; it can often be correctly determined which local minimum point corresponds to the 
interference of which pair of light waves. An example is the measurement of the thickness of 
a film for the purpose of quality control. In this case, plural sample films of the same kind 
are measured to determine whether the thickness of each film is appropriate, so that the error 
curves always take similar forms. Therefore, it is possible to empirically determine which 
local minimum point corresponds to the interference of which pair of light waves. In the case 
where some other method is available for measuring the thickness of each layer of a multiple 
layered film with certain accuracy, it is possible to estimate the positions of the local 
minimum points concerned from the measured values, taking into account the two-beam 
interference. 

The relation between the cycle interval T and the thickness d of the sample film for a 
local minimum point selected can be determined as follows. When, as is already known, a 
ray of measurement light enters a layer of refractive index n at an incident angle 8, the basic 
equation for calculating the thickness d of the layer is given by 

d [cm] = Ami {2 • (n 2 - sin 2 0) 1/2 (l/A 2 - Mi)} ... (7), 
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where Am is the order number of the peak of wavelength X 2 when counted from the peak of 
wavelength k\ toward shorter wavelengths. Normally, Am = 1. 

Using the cycle interval T [cm* 1 ] = \IX 2 - l^u the equation (7) can be rewritten as 

^=l/2-(« 2 -sin 2 0) 1/2 r ? 
and it is further rewritten as 

T = l/{2-(n 2 - sin 2 9) 1/2 d} = 1 X 10 7 /{2-(n 2 - sin 2 9) 1/2 d'}, 

where d is represented in cm, and d' in nm. 

When 9 = 0, the terms including sin 9 can be neglected, and the following equation is 

obtained. 

T=lX10 7 /(2-n-d')... (8) 
The equation (8) can be rewritten as follows to obtain the thickness of the layer. 

d'[nm] = lX10 7 /(2-n-T) ...(9) 

Substituting the cycle interval of the two-beam interference independently produced 
by each layer and the refractive index of each layer into equation (9) will give the thickness 
of each layer. 

An example of the above method for measuring the thickness of a film is described. 
The film in this example has two layers 21 and 22, as shown in Fig. 4. In this case, the film 
reflects three light waves R21, R22 and R23, so that there are three modes of interference: 
[R21, R22], [R22, R23] and [R21, R23]. The interference [R21, R22] is produced by the first 
layer 21, the interference [R22, R23] is produced by the second layer 22, and the interference 
[R21, R23] is produced by the combination of the two layers 21 and 22. Therefore, the 
thicknesses of the two layers of the film should be calculated by analyzing the first and 
second interferences [R21, R22] and [R22, R23]. 

Fig. 5 shows a graph in which the least square error calculated by equation (6) is 
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plotted against the cycle interval T. The least square error is calculated from the measured 
spectrum obtained by casting a monochromatic measurement light onto a two-layered film 
and scanning the wavelength of the measurement light. The graph shows two local minimum 
points at cycle intervals 12000cm" 1 and 14000cm" 1 . Though not shown in the drawing, there 
is another local minimum point at a cycle interval under 10000cm" 1 . Thus, the graph has 
three local minimum points. Each of the three local minimum points corresponds to each of 
the two-beam interferences [R21, R22] 3 [R22, R23] and [R21, R23]. 

The correspondence between the local minimum points and the two-beam 
interferences should be determined empirically, or by an appropriate method. An example of 
such methods is as follows. After the creation of the graph of Fig. 5, the thickness of the 
layer 21 is increased while that of the layer 22 is unchanged, and the measurement is 
similarly performed to create another graph. Comparison of the two graphs will show that 
one of the local minimum points has shifted in the direction of x-axis (i.e. wavenumber or 
cycle interval). This shift identifies the local minimum point that corresponds to the 
interference [R21, R22] produced by the layer 21. It is assumed here that the local minimum 
point at 12000cm" 1 corresponds to the interference [R21, R22] produced by the layer 21, and 
the local minimum point at 14000cm -1 corresponds to the interference [R22, R23] produced 
by the layer 22. From this assumption, it results that the third local minimum point at the 
cycle interval under 10000cm -1 corresponds to the interference between the reflected light 
waves R21 and R23 produced by the two layers 21 and 22. Therefore, the third local 
minimum point is not used in the thickness calculation. 

Letting di and d 2 denote the thicknesses of the first and second layers 21 and 22, and 
substituting the cycle intervals (12000cm -1 and 14000cm" 1 ) and the refractive indices ni and 
n 2 into equation (9), the thicknesses di and d 2 are given as follows. 

d,[nm]=lxl0 7 /(2 • m • 12000) 



16 

d 2 [nm]=lxl0 7 /(2 • n 2 ■ 14000) 

The method described in the Japanese Patent Application No. 2002-147107 also uses 
an error function defined similar to equation (6), in which thickness is used in place of cycle 
interval, and the error function contains reflective index. In contrast, by using cycle interval 
5 as a variable, the method according to the present invention uses no refractive index in not 
only equation (6) but also equations (4) and (5) to be substituted into equation (6). This 
means that the thickness of each layer can be measured even when the sample film consists 
of plural layers differing in refractive index and thickness. 

The foregoing discussion has assumed that the sample film is a multiple layered film. 

10 It should be noted that the method according to the present invention is applicable to the 
measurement of a single-layer film. In this case, the two-beam interference occurs at only 
one wavenumber, so that the least square error has only one minimum value. 

An apparatus embodying the present invention is described. Fig. 6 schematically 
shows the construction of the apparatus of this embodiment. 

15 The apparatus includes a light source 41, a spectroscopic section 42, an optical 

system 43, a photodetector 44, a spectrum-creating section 45 and a calculating section 46, 
where an interference spectrum created in the spectrum-creating section 45 is sent to the 
calculating section 46 and a thickness is calculated based on the interference spectrum 
according to the method described later. The spectrum-creating section 45 and the 

20 calculating section 46 are actually constructed by running appropriate programs on a 
personal computer. 

The basic operation of the apparatus is as follows. From the white light emitted 
from the light source 41, the spectroscopic section 42 extracts a ray of monochromatic 
measurement light having a single wavelength. The measurement light is cast on a sample 
25 S through the optical system 43. The sample S may be either a single-layer film or multiple 
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layered film. It is assumed here that the sample S has i layers. The sample S reflects the 
light on the front and rear surfaces, and also on the boundary surfaces in the case of 
multiple layered films. The reflected light waves travel again through optical system 43 
onto the photodetector 44. Receiving the light waves, the photodetector 44 generates 
5 electrical signals corresponding to the total intensity of the reflected light waves, and sends 
the signals to the spectrum-creating section 45. The reflected light waves interfere with 
each other, producing an interference light, as described later. When the wavelength of the 
measurement light is changed (or scanned), the output signal of the photodetector 44 forms 
a cyclic spectrum (or an interference spectrum). Based on the intensity signal from the 

1 0 photodetector 44, the spectrum-creating section 45 creates the interference spectrum with 
the wavenumber (or wavelength) as the abscissa and the intensity as the ordinate. The 
interference spectrum thus created is the measurement spectrum described above. 
Receiving the interference spectrum, the calculating section 46 calculates the thickness of 
each layer of the sample by a process according to the above-described method. 

1 5 The calculation performed by the calculating section 46 is described, using Fig. 7. 

In Step SI, the following information is externally specified as the calculation 
conditions: the upper limit 7a, the lower limit Tb and the interval A T of the cycle intervals to 
perform the calculation; and the refractive indices (&=1.../) of the layers of the sample. 
When the cycle intervals at which local minimum points will appear have been roughly 

20 estimated, the upper and lower limits Ta and Tb should be determined so that the range 
Ta-Tb covers the estimated cycle intervals, to reduce the calculation time. In the case of 
using a multiple layered sample, the calculation should cover the entirety of such local 
minimum points that correspond to the interferences each produced by each layer of the 
sample. Shorter intervals A T will provide higher calculation accuracy, but will require longer 

2 5 calculation time. Taking into account the trade-off between the accuracy and the calculation 
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time, the operator should determine the interval AT 

In the case of measuring a multiple layered film, Step SI may further include a 
sub-step of externally specifying an assumed range of cycle interval for each local minimum 
point corresponding to the interference produced by each layer. The assumed range is used in 
5 Step S8 described later. The assumed range of cycle interval for each layer may be otherwise 
specified in Step S8. 

In the case of repeatedly measuring the same kind of sample, it is possible to 
pre-register the information to be entered in Step SI so that the operator does not need to 
enter the information every time the measurement of the sample is repeated. 
10 When the calculation is started, the cycle interval T is set at Ta (Step S2), the 

constructed spectrum fa(x) or constructed vector Y' is searched for, and the least square error 
corresponding to the constructed spectrum fa(x) is obtained, as described above (Step S3). 

Next, it is determined whether the cycle interval T at the moment is equal to or 
greater than Tb (Step S4). When T is less than Tb, the cycle interval T is incremented by A T 
15 (Step S5), and the operation returns to Step S3. Steps from S3 to S5 are repeated until the 
cycle interval T becomes equal to or greater than Tb in Step S4. Thus, the range from Ta to 
Tb is scanned to search for constructed spectrums or constructed vectors at intervals AT, and 
the least square error is obtained for each of the constructed spectrums. 

The result of the calculation is represented by a curve showing the relation between 
2 0 the cycle interval and the least square error, as shown in Fig. 5 (Step S6). Then, this curve is 
searched for one or more local minimum points (Step S7). 

Within each assumed range of cycle interval specified for each layer in Step SI, the 
correspondence between the local minimum point obtained in Step S7 and the two-beam 
interference is examined (Step S8). That is, when an assumed range contains only one 
2 5 minimum point, it means that the minimum point corresponds to the two-beam interference 



19 

produced by the layer corresponding to the assumed range concerned. 

If a minimum point belongs to none of the assumed ranges, or if there are plural local 
minimum points within a single assumed range, it is impossible to automatically determine 
which two-beam interference corresponds to the local minimum point. In such a case, based 
5 on the curve showing the relation between the cycle interval and the least square error, the 
operator should determine, empirically or by other methods, which local minimum point 
corresponds to which two-beam interference, and enter information indicating the 
correspondence into the calculating section 46. Alternatively, it is possible to remove the 
step of specifying the assumed ranges from Step SI and instead provide Step S8 with a step 

10 of letting the operator manually determine the correspondence between all the necessary 
local minimum points and the two-beam interferences, based on the curve showing the 
relation between the cycle interval and the least square error. 

To calculate the thickness of each of the first to /th layers, the calculation by equation 
(9) is performed, using the cycle interval corresponding to the local minimum point 

1 5 identified in Step S8 and the refractive indices entered in Step SI (Step S9). 

The above embodiments are mere examples of the present invention and may be 
modified in various ways. For example, the base spectrums used for representing the 
constructed spectrum may be differently defined. When the film is very thin, the waveform 
of the interference spectrum tends to have peaks or troughs whose cycle interval cannot be 

2 0 clearly determined. In such a case, the cosine terms and/or (1/x) terms may be removed from 
the base spectrums to obtain a better approximation result by easier calculation processes. 
Furthermore, the approximation accuracy can be improved by appropriately changing, 
adding or removing a term or terms of equation (3'). 

In the above embodiments, the interference produced by two or more layers 

2 5 combined is not taken into account because the thickness of each layer can be calculated 
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from only the interference produced by the layer concerned. However, when the interference 
produced by each layer cannot be separately extracted because of plural local minimum 
points located too closely, the thickness can be calculated by using the interference produced 
by a combination of two or more layers. For example, in Fig. 1, when the local minimum 
point corresponding to the interference [Rll, R12] is located too close to another local 
minimum point to calculate the thickness dj of the first layer from the interference, the 
thickness dj can be calculated by 

T, 2 =lxl0 7 /(2-(n 1 di+n 2 d 2 ), 
where d 2 is the thickness of the second layer calculated from the interference [R12, R13], 
and T12 is the cycle interval at which the local minimum point corresponding to the 
interference [Rl 1, R13] appears. 

As described above, the method and the apparatus according to the present invention 
greatly improves the thickness measurement accuracy by taking into account various factors 
appearing in the measured spectrum, such as the offset, drift, phase shift, and amplitude 
change of the cyclic waveform. Also, the method and the apparatus according to the present 
invention makes it possible to measure the thickness of each layer of multiple layered films 
even when the layers differ from each other in thickness and refractive index. 



